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1. Introduction 

The concepts underlying tlie Luttinger model (LM) [HEllSlll! (for reviews cf. [SI El 13 IE]) 
find increasing application in the theory of quasi one-dimensional ultracold quantum 
gases. Using the harmonic fluid model of |9^, properties of a one-dimensional trapped 
Bose gas were investigated in pLOJ. Related concepts were also used in the theory of 
interacting bosons in one- dimensional optical lattices fUE]. In ^3], mass-composition 
separation as analogue of spin-charge separation in the LM was studied for a two- 
component neutral Fermi gas in a one-dimensional harmonic trap (cf. also Jl]). A 
detailed description of the harmonic fluid approach to one-dimensional quantum gases 
is given in ^3]. 

The present authors formulated a theory of interacting one-dimensional fermions in 
a harmonic atom trap [T^ IT7j based on bosonization in analogy to the LM. The theory 
was initiated by the prospect to create experimentally a highly degenerate Fermi gas of 
neutral atoms in a quasi one-dimensional harmonic trap. A possible route to achieve this 
aim is microtrap technology f cf . |18t IT^ 1^ 1^ I22j ) . This will allow to study non- Fermi 
liquid behavior of interacting one-dimensional Fermi gases without the complications 
due to contacts and impurities. 

The model in E|; which can be termed Tomonaga-Luttinger model with 
harmonic confinement, lacks so far a phase operator formulation known for many models 
of one- dimensional interacting fermions. Probably best known for such a phase operator 
formulation is the Tomonaga-Luttinger model with periodic boundary conditions. A 
phase formulation for the corresponding model with open boundary conditions (OBC), 
which belongs to another universality class is found in [21] and in particular in 
. Related phase theories also exist for the case of interacting one-dimensional bosons 

pea in. 

In general, a phase operator approach to one-dimensional fermions comprises several 
features: 

i) Extending the linearized dispersion near the Fermi energy with respect to fermionic 
creation and annihilation operators c^^ and Cka to all values of k and adding an 
anomalous vacuum of filled negative energy states, Kronig's identity [22] can be 
used to transform the free Hamiltonian into a free bosonic form involving phonon- 
like density fluctuation operators 6+^ and bg^^. The index cr = ±1 accounts for two 
spin directions or two different components in the Fermi gas. 

ii) The pair interactions between fermions must be describable in terms of bilinear 
products of density fluctuation operators. 

iii) In the spinful case, an additional transformation to charge and spin degrees of 
freedom (or mass and composition degrees in the neutral case), denoted by z/ = ±1, 
is needed on the way to a diagonal Hamiltonian: 

dqu = A=^(^^ bg^, i^ = ±l. (1) 
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iv) In the absence of backscattering between different components, the Hamiltonian 
then separates into decoupled spin and charge (or mass and composition) parts 
which can be diagonalized by a Bogoliubov transformation. The final Hamiltonian 
H for the excitations of the one- dimensional Fermi sea typically has the structure 
H = J2qu '^qu \q\ fqufqt^ terms of free phonon operators /"*" and / and renormalized 
velocities Vqy. Zero modes |lj must be added to make the correspondence with the 
original fermionic Hamiltonian complete. 

v) If the dependence of the velocities on g, i.e., that of the underlying coupling 
functions is neglected, hermitian phase fields ^u{z) and dual fields Qu{z) can be 
defined in terms of the ^-operators (or else in terms of the /-operators) which 
contain zero modes and bring the Hamiltonian into the canonical form 



^ = E y jd^ [^K, {dM^)f + ^ {dMz)f] , (2) 

where are the central coupling constants. Backscattering between the two 
components Vj_ oc f dz cos \/S^^i{z) destroys the simple quadratic form of the 
spin part [771I2H]- 

vi) An important step in any phase theory is the representation of fermion field 
operators in terms of the phase operators in order to calculate correlation functions 
easily, a procedure usually called "bosonization" . The method goes back to [22] 
and [SO! and is now well understood |SI1 ll| • For open boundaries, it is described in 
the works EH and 



Here, we consider interacting one-dimensional fermions in a harmonic trap. Then 
the single particle wave functions are not simple combinations of plane waves. We, 
therefore, use a general bosonization prescription described in |22] which involves 
auxiliary Fermi fields "ipaa- related to the original operators Qo- (note that for the harmonic 
oscillator q is replaced by a discrete index I = 2,— 1,0,1,2,...) by 



oo 



V'..(m) = E e*''^ = iJaAu + 271). (3) 

/=— oo 

These operators are defined on the unit circle as will be most other operators below. 
Bosonization is provided by 



in terms of the bosonization phases 

M = E ^ e^-(«+^./2) (5) 
^ 'm 



m=l 

The quantity t] in the non-normal ordered expression is a positive infinitesimal. 
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There is a nontrivial relation between the bosonization phase + 0o-('u) and the 

physical phase operators ^u{u) and 6y(u) given below. The Klein operator Ua- reduces 
the fermion number by one. An explicit construction of k^r = — ilnf/^ is given in |32j . 
The zero modes in the physical phases and Bj^ are related to uN„ and ko- by the 
transformation equation (^. 

In the present case of the harmonic trap, the relation between the auxiliary variable 
u in the correlation functions and the spatial position z inside the trap turns out to be 
u uo{z) = arcsin(2;/Li?) — 7r/2, where 2Lf is the quasi-classical extension of the Fermi 
sea. 

The paper is organized as follows: Section 2 develops the phase theory for a one- 
component Fermi gas according to the above prescription. Section 3 extends the theory 
to two components. In Section 4, we calculate the scaling exponents contained in the 
one-particle correlation functions. The Appendix gives all details of their evaluation 
using the WKB representation developed in 



2. Phase theory for one component 

We start with the simpler case of one fermionic component, e.g., spin polarized fermionic 
atoms in the case of ultracold atoms. The one-dimensional equivalent of s-wave 
scattering, i.e., a contact interaction is forbidden by the exclusion principle and the 
interaction in such a Fermi gas is usually weak. The model may, nevertheless, have 
potential application to ultracold Fermi gases in a harmonic trap provided Feshbach 
resonances enhance the interaction as recently found in the three-dimensional case |35j . 

Quasi one- dimensionality can be realized when the particle number N is less than 
the ratio of transverse to longitudinal trap frequency uJt/uJi and the three-dimensional 
interaction energy per fermion is less than the transverse excitation energy in the highly 
elongated cylindrical trap jSH]- 

However, it is not enough to reach degeneracy, ksT <^ ep- In order to observe zero 
temperature behavior, the much stronger condition ksT -C huji is also required. 

The model considered in ^\ is described by the bosonic excitation Hamiltonian 

H =^huji ^ m + + ^Y1 [^cim) - Vb{m)] V2m 4™ + dtm (6) 

+ 1 E ^a{m) m [d^di + d^djj - ^ ^ [V,{m) - m (d^ + d^^ 

In the case of one component, one does not have to discriminate between d and b 
operators. We treat the one-particle operators from Vb and Vc on the same footing to 
avoid problems with spurious self-energies when Vb is retained. 

The interaction terms in equation (jHl) follow uniquely by retaining those parts in 
the fermionic pair interaction operator 
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V = V{m,p;q,n) {c^Cg){c+Cn), (7) 

^ mnpq 

which are expressible in terms of density fluctuation operators p(m) = J^pCp+m^p (in 
the basis of harmonic oscillator states ipniz)), i.e., 

V{m,p; q, n) Va{\q - m\) 5m-q,n-p + - m|) 6g-rn,n-p + - P\) Sm+g,n+p- (8) 

In jnZ], it is shown that the retained matrix elements are dominant in the one- 
dimensional many-fermion system in the harmonic trap. This is related to approximate 
momentum conservation during collisions in the trap. 

The interaction matrix elements Va{m), Vh{m), and Vc{m) correspond to the 
Luttinger model coupling functions (74 (p), (72 (p), and gi{p), respectively. Va and H 
describe forward scattering and Vc describes 2k p (backward) scattering. The notion of 
backscattering in the one-branch-one- component system is to be interpreted as follows: 
The study of the matrix elements showed that the right (left) running part of a one- 
particle wave function interferes with the left (right) running part of the wave function 
of another fermion exchanging a momentum of about 2k p. 

The drastic reduction of interaction matrix elements necessary to obtain equation 
(jni) from the full Hamiltonian for a given one-dimensional (effective) pair interaction 
between fermions may affect physical details. However, the model will still provide us 
with the correct critical exponents in terms of the central coupling constant K. This 
has been explicitly shown in for OBC. It must be pointed out that the model 
accounts for the basic interaction mechanisms in one dimension, namely forward and 
backward scattering. It may be more practical to consider K as a parameter taken from 
experiment, as is usually done in applications of the Luttinger model. 

2.1. Phase operators 

The second contribution on the r.h.s. of equation (jH]) represents a one-particle operator 
Vi. This operator is neglected in the model with open boundary conditions as pointed 
out in inSj. The one-particle operator does not alter the critical behavior, but has 
quantitative effects on other properties. Vi is exactly taken into account in the present 
model. For more details, we refer to [T7] . 

The Bogoliubov transformation to diagonalize the Hamiltonian in terms of 
operators / and /+ is given by = fm cosh Cm - sinh + r]^ exp(-Cm) , where the 
inhomogeneous part r]2m = —[Vdrn) — Vb{m)] exp(— ^2™)/ (2e2m"\/2m) originates from Vi 
and differs from zero only for even indices. 

The central dimensionless coupling constants Km, transformation parameters ^m, 
and the renormalized level spacings for the model, equation are given by J7j 
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-Kr, 



\ huje + Va{m) + {Vb{m) - V,{m)) 



2n 



(9) 



respectively. 

In the phase formulation, the dependence of K and e on m must be suppressed. 
Note that Va{m) Va implies that there are no interaction effects due to this 
matrix element in the one component theory and Va should be omitted, i.e., we use 
e = J{huJey -{Vb-V,y and 



K' = 2tiK = 27re'^^ = 2n 



hue - jVb - Vc) 
\ hue + {Vb - Vc) ^ 



hujf 



2K 



(10) 



where H is usually small in comparison to Vc, even in the marginally long-ranged case 
of the dipole-dipole interaction j37j. 

In some physical quantities, e.g., the one-particle correlation function, neglect of the 
dependence on m in i^' (which for large m must approach unity |3]) leads to divergences. 
We take care of this by writing 



1 + {K -1) exp( 



-mr 



'111 



with a small positive number r <^ 1, estimated as r R/Lp oc l/y/N where R is the 
range of the interaction. This small number should be discriminated from the positive 
infinitesimal t] appearing in equation For m < 1/r 3> 1, there is practically no 
dependence on m. 

We begin with rewriting the one-particle operator Vi. This already leads to 
the structure of the main phase operator. Assuming the same exponential decay as 
in equation (jllj) also for interaction matrix elements, i.e., V^(m) — Vb{m) = {Vc — 



Vb) exp( 



-mr 



Vi can be expressed as 



where 



V = -{V-Vb) 



du 



-r+2iu 



-r—2iu 



+ 



1 - e-^+2m l_ ^ 



du<Podd{u), (12) 



<Podd{u) = - (0(m) + 0+(m) 



-u 



\-u)) = J2 \l— sm{nu) [dn + d:^] . (13) 

n=l 



n 



The same phase operator appears in the particle density operator when it is decomposed 
into a slowly varying part and a part describing Friedel oscillations 
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with kq = {K'^ — l)/8. The main phase operator corresponds to that in 



°° Tif p—mrj 

$(n) =Nu + <Podd{u) + h{u) =Nu+Y. \ — sin(mn) (/„ + /+). (15) 



m=l ' ^ 



Note, that in [23] rescaled phases are used and the function h{u) is absent due to the 
neglect of one-particle operators. The dual phase operator is 



= ^ E y cos(nM) [dn - dl) = —Y.\l cos(nu) (/„ - /+j . (16) 



We prefer not to combine the zero mode operator k with the dual phase field 0. 
Instead, we retain the unitary operator U = exp{ik) in order to avoid the mathematical 
problems pointed out in The relevant commutator is [exp{ik),N] = exp{ik). The 
phase operator corresponding to the momentum density is 



• oo 



n(u) = duQ{u) = V Vn e-"'' sin(nu) (d„ - rf+) (17) 

2^ n=l 



n=l 



Using E^=i (e*™("+^''/2) + ^-im{u-ir,/2)^ _^ 2TT62n{u) - 1 Icads to the commutator 

<I(m), Il{v)] = ^ -v)- + v)) , (1^ 



which is also implicit in the work [25]. This is not a canonical commutator. Canonical 
commutation relations cannot hold for operators with fixed parity. 

2.2. Phase Hamiltonian 

Inserting equation (fTH|l and equation (fT7]l into 



H = - 



2 J-^ 



du 



K'Ii\u) + — [du^u 



duH, 



(19) 



and neglecting an irrelevant contribution from b'^{u), reproduces equation ^ in addition 
to the zero mode contribution Eo[N) 
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EoW = ^N' = j^^N'. (20) 

In the latter relation, corrections 0{N) are neglected. Remarkably, canonical equations 
still hold in spite of equation (fT8|) : Consider the Heisenberg phase field t) (we set 
^ = 1) 

idMu, t) = [I>(m, t), H{t)] = ie K' fl{u, t) = i—. . (21) 

dll{u,t) 

Thus dt^{u,t) = e K' fl{u,t) remains correct. 
2.3. Bosonization 

A simple expression of the individual operators (f) and (f)'^ in equation in terms of 
the new phase operators $ and 6 is not available. However, their sum which appears 
in the non- normal ordered form of equation (j3)), is represented by 

0(m) + 0+(m) = ^oddiu) - 27r 0(m) = $(«) - b{u) - Nu - 2n 0(m). (22) 
This can be utilized to calculate the auxiliary correlation function 

with the chemical potential 



/i^ ^ Eo{N) - Eo{N - 1) = [n - . (24) 



Analogously 



Derivation and evaluation of these expressions is described in the Appendix. 
3. Phase Theory for two components 

In terms of mass and composition operators, i.e., after the transformation equation (^, 
the excitation Hamiltonian for forward scattering is given by 



dmu ~l~ d^jjd^J^ (26) 

mu 

+ W HVh (m) + z.H±(m)] {d-^Z + dl,] ■ 



2 m,v 
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It must be supplemented by backward scattering and one-particle contributions. Noting 
that the one-particle operator | Y.m>o,a V2m {Vc\\{m) — Vm{m))[b2ma + ^27710-]) which 
originates from rearranging the two-particle backscattering operator and from a spurious 
self-energy in V^, is transformed by means of J2a = 1 + v = 25y^i and 



the final form of the additional operator is 



(27) 



VaM = - 9 E rnV,\\{m) + d^} + ^ - + ^mi) ■ (28) 

m>0,i^ m>0 

So far, we did not mention Vcx, i.e., backscattering between the two components, 
but we note that Vci_ does not generate a one-particle potential. 

In omitting the dependence on m, Va\\ drops out again. Then the renormalized level 
spacings and central coupling constants are given by 



(29) 



and 



{hue + lyVa^) - + - Kii] _ Kl 



\\ {hwt + uVai) + + i^H± - Kii] 27r ' 
respectively. The total one-particle operator from T4 and can be rewritten as 



(30) 



-r+2m 



1 - e-''+2*« 1 - e- 



r—liu 



du^u,oddiu), (31) 



with 



0^.,odd(^i) = E ^ e""^/^ sin(nii) {dn^ + d+J . 

n=l V'^ 



(32) 



The corresponding phase shifts are the c-number functions (real and odd) 



Oi/(ti) = ? — ^ (1 + Jy) In 



— r+2iM ^ 



1 -e- 



-r— 2m 



mi 5„ 1 In 



1-e 



-r+2m ^ 



1-e- 



-1 — 2m ; ' 



(33) 



87re^\/2 

with Ki = (A'^ — l)/8. Because of 6_i = only the mass part is affected by the 
one-particle operator. The phase operators are 



M'^) = <t>u,oM{y) + Ki'^) +uK-Y.\^ 6--"/^ sin(nii) (/„, + /+ ) + uK, (34) 



n=l 



n 
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©^-l") = ^Y.^ COs(nM) (dnu -dnu) =^Y.\J COs(nM) {fnu - fnu)- (35) 

The dual phase fields Qy generate associated momentum densities via 

Il,{u) = — e,{u) = -— J2 ^ne~^^ sm{nu) - c/+ ) . (36) 
Finally, the phase Hamiltonian is 

^=E| jy^[K'Mu)+^{duMu)yy (37) 

The zero mode in gives the ground state energy operator 



Eou = j^\6,,iY.NaN^' + 6.,-i{Y.aN^] |. (3^ 

Again, the equations of motion for the Heisenberg fields ^y{u,t) are dt^u{u,t) = 
Bosonization according to equation (jD) is provided by 

0^(m) + 0+(tx) = ^ ^ - K{u) - uK - 271 Q,{u)). (39) 

Backscattering between different components is described by the operator Vc^. The 
composition part Hy=_i acquires the form of a sine-Gordon Hamiltonian on the unit 
circle. Considering the results for the usual sine-Gordon system, one can speculate that 
a sufficiently repulsive interaction renders Vc± irrelevant at low energies. Restricting the 
treatment to long range interactions as discussed in avoids this problem. 

4. Critical exponents 

This section compiles the results for critical exponents extracted from the one-particle 
correlation functions which are calculated in the Appendix. Subsections are devoted to 
the static two-point correlation function, the local dynamic correlation function, and to 
the local spectral density of states, respectively. 
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C{z^, Z2) = C{z,, t = 0; Z2) = {^+{z,) tij{z2))o. (40) 

We have to evaluate equation ()A.11|) using the explicit result ()A.17|) for t = 0. To this 
order, we introduce the abbreviations Z,^ = (1 — zl/LpY^'^; u = 1,2, Lp = t\/2N — 1 
[i is the oscillator length) to find ioY Lp — \zv\ > Lpr^ and zi 7^ Z2 (the case zi = Z2 
was considered in |3^ ) 



\Vi±2uoiz,)+tr)\' ^ -—, \Viuoiz,)±Uoiz2)+tr)\' ^ , ,^ , -j-. (41) 

Z [L ± ZiiZi2 — Z1Z2/ iJp) 

This gives 

sm{kF{zi - Z2) + A_) 



C{zi, Z2) 



^27o2ao + |-70 



N-1 cos{kp{zi + Z2) + 



(42) 



(1 + Z1Z2 - z^z2/Liy'"^-2{i - Z1Z2 - Z1Z2/L 



{Z1Z2 



i"o- 



with the (irrelevant) phase shift A± = | (arcsin(2;i/Li?) ± arcsin(z2/-^^F)) — b{uo{zi)) =f 
b{uo{z2)) and notations from the Appendix, e.g., z^, = z{z^) = z^Z{z^)/2 + 
{Lp/2) aTcsm{zu/Lp). For zi = z > i and Z2 = 0, and well inside the trap {\z\ ^ Lp), 
we find 

^centeAz) - ^z/Lp)W + ^ 2tiLp {z/Lpf^'^- ^ ' 

Because of \z\/ Lp ^ 1 and 270 + 1 = 2ao + K, the second term is not important. Thus 

j^^'Yo siri Ivjpz 

Ccenter(^) = ^ (^/L^)27o+i' (^^^ 
gives the decay exponent (twice the bulk scaling dimension of ip) 



= 270 + 1 = ^ 



2A. (45) 



This result was obtained numerically for the present model in the second paper of |17j . 
It coincides not unexpectedly with the OBC result |211I23ES] in the bulk limit and also 
with that for the homogeneous Luttinger model. 

We can extract from equation the boundary exponent A^. To this order, the 
thermodynamic limit must be performed explicitly: As in jHS] this is done by making 
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the trap shallower and shallower: Ui oc — » 0. Furthermore, the coupling constant K 
and the Fermi wave number kp must be held fixed. For a specific interaction model, this 
limit was already performed in Furthermore, besides r] also r oc vanishes 
and one should consider renormalized fields as described in I^U] . 

Consider then the boundary situation Z2 = Lp{l—ri) while < zi Lp, such that 
Z2 —>■ Lp that is ^2 — > in the denominator of C{zi,Z2). The distance Z2 — z\ then is 
large as required for Aj^. Equation gives 



Lp{\ - f])) OC ^ r iVvr''^' oc If^^li ^. (46) 



Hence 



1 3 1 

2Ax = - ao + 7o + ^ = g 



(47) 



This boundary exponent is different from the corresponding result for OBC in 
4-2. Local dynamic correlation function 

Following we consider only the part varying slowly in space, i.e., we neglect the 

rapidly oscillating Friedel part given by the last two terms in equation ()A.11|) . We then 
obtain 

C''^{z,t;z) = ^^^^^^^ {i^nuo{z),t)Muo{z)))o. (48) 
According to equation ()A.17|) . 



{ij^{uo,t)iJaiuo))o = e^^'r^-^^ \v{-et + irf^''V{-et + ir,)\ (49) 

f V{2uq -et + ir) V{-2uo - et + ir) 1 "° 
I \V{2uo + tr)\^ J 

holds with uq = aicsinlz/Lp) — tt/2. The renormalized level spacing e vanishes in the 
thermodynamic limit, but the chemical potential becomes a constant independent of 
particle number: fi^ = e {N — 1/2)/ K ^ 2h/{K'^ + 1) {ujeN) = const.. Vanishing rj and 
r imply 

^<-'*' = -2ii^(575y ^'2"»-'*'^<-2"°-'*' = 4M575pT37iy- (50) 

In order to collect all possible powers of t, we specify z to be so near the classical 
boundaries that cos^(et/2) — z^ / L"]^ oc t^ holds. Then we obtain: 
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C'' oc V{-et) 
with the boundary exponent 



27o+2ao+l 



1 



OC 



t 



2A|| ' 



(51) 



2A,| = 270 + 2ao + 1 



1 



(52) 



This agrees with the one-component results in ESI IBHI for OBC. 

However, for z away from the boundaries, only D(— et)^'''"^^ contributes and the 
bulk scaling exponent is reobtained 



270 + 1 



1 
2 



K + 



1 

K 



2A. 



(53) 



Note that the scaling relation 2A± = A + Ay ^41, is not fulfilled with respect to the 
variables z and t. However, boundary conformal field theory |12] applies to the auxiliary 
model which uses the variables er {it = r) and u: The infinite strip w = er — iu with 
— TT < u < can be mapped conformally onto the complete complex plane and the 
Euclidean Lagrange density derived from equation (fT!I|) is locally rotation invariant. 
Consequently, the scaling relation holds irrespective of the boundary condition. This 
can be checked by direct calculation: If we calculate (V'^(m ~ —TT/2)ipa{v ~ 0,7r))o 
instead of the physical correlation function C{zi ^ 0, 2:2 ~ Lp, —Lp) in order to get the 
boundary exponent A^,^'' for the auxiliary fields we find 



{iptiu) '^a{v))o OC 



sm 



2A 



(a) ) 



(54) 



with the new value 



2A 



(a) 



1 

4 



K + 



K 



(55) 



The critical exponents A and Ay remain unchanged so that the scaling relation is 
fulfilled. 

The subsequent transformation to the physical plane = z + ier {z = z/Lp) which 
must ensure z = cosm is, however, not conformal. The breakdown of local conformal 
invariance with respect to the scaled coordinates z and er near the boundaries becomes 
obvious by inspection of the corresponding Euclidean Lagrange density 



(56) 



with 
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We must evaluate the anti-commutator omitting the Friedel part 
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NF 

A{t,z) = ( ^+{z,0),^{z,t) ^ ) . (57) 







Using equation ()A.17|) and equation ()A.19|) of the Appendix, we find 



w X / , iV(-2un-et + ir)V(2un-et + ir)V\ 
A{t, z) = ( e-^^-+^*I?(-et + irf^-V{-tt + ir^) \ p^2u,)\^ ^ 



+e--- Viet + .r)^- V^et + | Hz^+ + -)^(2-o + et + ^r) \ ^ \ r^- 



|r'(2Mo)P J / -^LpZ' 

In the thermodynamic hmit, r — 0, /ia? — /iat+i — /i, and measuring the frequency uj 
from the chemical potential, u — ^ u, we obtain 

N{cu, z) = — dt A{t, z) oc y (59) 



oo 



^ -i(70+l/2)7r+i(7,)+l/2+ao)et 

X / dtcos(ujt) { — „ , -I , — — - — -- — — ^ , „ + CO. 

lo ^ '[ sm^^°+\te/2) [cos^{te/2) - z^/Llf' 

We further proceed as in The regularization cos(co't) [cos(co't) — 1] is adopted 

and the integral f^ dt[cos{t) - l]/t^ = 7i/{2r{k) cos(7rA;/2)), 1 < A; < 3 is used. Well 
inside the trap, we find 

' ' (60) 



Ar(.J,|<<i,)«____y 
The bulk exponent of the local spectral density thus is 



1 

"bulk = 27o 



< 70 < 1, (61) 



2 

and agrees with the corresponding result in flEi ES] for OBC 

The boundary exponent is obtained as follows: Near a boundary, \z\ ~ Lp, we 
obtain for iZ^ < e^t^ < when [cos2(te/2) - z'^/Llf 6*^02-2^0 (ig)2ao^ 

iV(^, N ^ OC ^ ' \^ ^ZJ -\ . (62) 

vrLiT er(27o + 2q;o + 1) Ve/ 

The boundary critical exponent thus is 

1 , . 

"boundary = 270 + 2ao = — - 1, < + 7o < 1- (63) 

Again, this value agrees with the one-component version extracted from [211 123 IHH] for 
OBC. 
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4-4- Two components 

It is straightforward to extend these results to the case of two components: The 
correlation function factorizes in the mass and composition representation. However, 
the corresponding exponents in equation ()A.12|) have half the previous value due to the 
transformation equation This leads to 



2Aii 



1 
2 



1 1 

Ki ^ K^i 



(64) 



2Ai_ = — 
16 



Ki K-i 



For the spectral functions, we assume equal particle numbers of the two components. 
Then the spectral function calculated above (and amended by the correction just 
described) refers to each component giving 



"bulk 



1 

4 



'^boundary 



1 
2 



1 1 

^1 ^ H- 



(65) 



These values agree fully with those found in [211 12^1 IBS]; except for the value of A^. 
Again, this results from the violation of local Lorentzian invariance with respect to 
spatial position and time in the present model. 



5. Summary 

We gave a consistent phase operator formulation for a model of interacting one- 
dimensional fermions in a harmonic trap. The model is the analogue of the Tomonaga- 
Luttinger model with open boundary conditions (OBC). Suitable phase operator were 
identified by expressing the one-particle operators as linear odd combinations of the 
basic bosonization operators in agreement with earlier results on the decomposition of 
the particle density operator. The total Hamiltonian became a quadratic form in the 
phase operator $ and the operator fl corresponding to the momentum density in spite 
of the presence of one-particle operators and the harmonic trap potential in the original 
fermionic Hamiltonian. These results were extended to the case of two components. 
The phase formulation bears many similarities to that for OBC. However, distinctive 
differences are also present. The spatial position is non-linearly related to the variable 
in the phase fields. This spoils local Lorentzian invariance usually satisfied by related 
phase theories. In addition, form factors appear which take account of the harmonic trap 
topology. Exact results for the one-particle correlation functions were derived and used 
to extract bulk and boundary critical exponents including those for the local spectral 
density. The values of the critical exponents coincide with those for the OBC except for 
the boundary scaling dimension of ip. 
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Appendix: Calculation of one-particle correlation functions 

5.1. Bosonization and auxiliary correlation function 

We use the representation of auxiliary Fermi operators |32] 

^+(u) = 6+(m) e-*<^^(") e-*^(") = — t/+ e-<'^°''''(")+^"~2"®(")) . (A.l) 
The Klein operator 0{u) = exp{iNu) U acts on the fermion number operator N. U 



commutes with all bosonic operators 







etc. and changes the 



fermion number according to f{N)U'^ = U^f{N + 1). The contribution Nu to the 
full phase operator in equation (fT^ leads via equation (fT^ to the ground state energy 
equation (pn|) and consequently to the chemical potential equation (j2H). We consider 
the one-particle correlation function of auxiliary Fermi operators 



Ca{u,t;v)= N{'lpti'^^t)'^a{v,0))N, 

evaluated with respect to an N-fermion state. The bosonized form is 



(A.2) 



1 

V 



Ca{u,t;v) = - N{d+{u,t) e-^('^+(-'*)+^(-'*))e'(<^^(")+^(''»0(i;)) 



N- 



The time evolution operator is governed by the separated Hamiltonian 



(A.3) 



H = H + EoiN) =J2memf;l,fm + EoiN), 



(A.4) 



made up of the ground state energy of the N-fermion system plus the Hamiltonian of 
collective excitations. We thus get 



In order to determine U'^{t), we set ^ = 1 and f{N,t) = exp iEo{N)t 



and find 



U+{t) = f{N, t) U+ f*{N, t) = U+ f{N + i, t) f*{N, t) 
= f/+ exp [i (^Eo{N + i) - Eo{N)^ t] = f/+ exp [i/i^+it 



(A.6) 
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Thus the time dependence of JJ^ is provided by the chemical potentiah U^(t) = 
exp{ifij^t) tj^ . The final result for Ca is: 

Ca{u t v) = 1 e"*^^^"^'"'"^^-'^^^^* ^g-*(^odd(«.*)-27r0(u,t)) gi(0odd(f )-27rG{t;)) ^ 

The remaining expectation value is purely bosonic and can also be thermal. 
5.2. Calculation of one-particle correlation function 

The next step is to actually calculate the physical one-particle correlation function 

e^"^{i;+{u,t)Mv)) (A.8) 



This expression contains sums over harmonic oscillator wave functions ipmiz). In 
a WKB method was proposed to deal with such sums: For large N, the (singular) 
expansion 

was found. Here, the following abbreviations are used 

/Z\7l_l l Z 

Uo{z) = arcsin ( — I — — , z = -zZ[z) + -Lp arcsin -— . (A. 10) 

\Lp / 2 2 2 Lp 

The spatial coordinates Zi and Z2 are thus restricted to the classical region {—Lp,Lp). 
The correlation function takes on the much simpler form 

an / 4 \'/^ 

x|e*'=^(^(^^)-^'(^^))C',(-Mo(zi),t; -Mo(^2)) + e-^'=^(^(^^)-^'(^^»C,(Mo(^^ 



Phase theory for confined ID fermions 



18 



Note the similarity with the representation of the one-particle correlation function for 
OBC in • However, the argument in the operators is non-linearly related to the 

spatial positions due to the harmonic trap. In order to apply the Wick theorem 

(e^e^) = exp ((i^) + ^{A') + 1(5^)) , (A.12) 

the operators A and B must be homogeneous linear combinations of / and f'^. We call 
the corresponding part in equation (fTK|l <!>/, hence 0odd = $/ — &, and get 



(A.13) 



Thus A = -i^f{u, t) + 2TiiQ{u, t) = A{u, t), B = i^f{v) - 27iiQ{v) = -A{v, 0). 

We begin with the case of zero temperature when only {fmfn)o = Sm,n survives. 
Then simple but lengthy calculations using equation (fTT|) and l/K^ — * 1 + (1 — 
K) exp{—mr)/K and the coupling constants 



K 



- K 



7o 



give 



exp ( ^{A\u))o ) = v^r^" {V{2u + ir)V{-2u + ir)} 



(A.14) 



(A.15) 



with Vi^s) = 1/(1 — exp(is)). The correlation function {A{u,t)B{v))o 
— {A{u,t)A{v,0))o turns out to be 



{A{u, t)B{v))o = lnV{u -v-et + iri) 



(A.16) 



+7o {lnT>{u — V — et — ir) + \nV{v — u — et + ir)} 



+ao[\nV{u + V — et + ir) + \nV{—u — v — et + ir)]. 
The final zero temperature result is 



Ca{u, t- v) = t)iJa{v))o = e-^(^-i)(«-^)+^'^^'*+*^(«)-^^(^) r^T" (A.17) 

[D(u — V — et + ir) T)(v — u — et + ir)]^° j T)(u + v — et + ir) T)(—u — u — et + ir) 1 "° 
^ 1 - e-'?+^(«-^-^*) \ \V{2u + ir)V{2v + ir)\ J 

Note that in contrast to fl^ two small parameters appear in equation ()A.17|) . the positive 
infinitesimal rj and r which is finite in the finite system. In addition, the one-particle 
operators in the Hamiltonian lead to additional phase factors which can be expressed 
^g. ^i(b{u)-biv)) ^ {x>{2u + ir) V{-2v + ir)/{V{-2u + ir) V{2v + ir))}''\ 
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The main difference of equation ()A.17|) to the formula (25) in j^HI is the appearance 
of the auxiliary variable u (or Uo{z)) in place of the position. Otherwise, our result ()A.17|1 
is in direct correspondence to that formula. As a consequence, the finite temperature 
calculation of |^ can be adopted to obtain the result for the canonical ensemble by 
making the substitution 



V(s) 



1 



1 - 



1 + 



sin(s/2) 
sinh{k(3e/2)^ 



(A.18) 



with /3 ^ = ksT. We will not pursue the finite temperature case. In the same way, we 
find 



{^a{u, t)^tiv))o = e^^(«-'')-^^'^+i*-*^(")+*''('') r^^o (A.19) 
['D{u — V — et + ir) T){v — u — et + ir)y° [ T)(u + v — et + ir) T>{—u — v — et + ir)\"° 



^ g—T]+i{u—v — et) 



\V{2u + ir)V{2v + ir)\ 



Appendix B: Backscattering between the two components 

In the approximation Vc^(|m|) Vc±, the associated two-particle operator in 
second quantization is 



Using 



X /""' dv ^ ^ ^ 



(v) = e-*(^--i)'' ij+ —(,-i4>t{v)-i4>a{v) 



(B.l) 



(B.2) 



together with equation (jH), and utilizing \tj(j,N^i] = 5(jytJ„, we obtain the following 
phase representation of V± 



V± = Yl J c?Mexp[2m(A/'^ - A^„^)] e^^[2i{ip„^odd{u) - (p-^^oddiu))], (B.3) 



with new phase operators, 

oo 



2m 



sin(m-u) {dmy + 



u m=l 

Because of ipa,odd{u) — 0-a,odd{u) = (Tv^(<l'_i — uNu=-i) the simpler form 



2T^rf J -TV 



du cos 



u] 



(B.4) 



(B.5) 



is obtained. As expected, only the composition part is affected by backscattering 
between components. 

Finally, it is noted that that normal ordered form of equation ()B.5|) is 

V8I'_i(m) 



du cos 



TT 27r ■ 1 + e-^*? - e-'' cos(2m) 



(B.6) 
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